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Abstract 

We study the distribution of optimal path lengths in random graphs with random weights as- 
sociated with each link ("disorder"). With each link i we associate a weight n = exp(arj) where 
rj is a random number taken from a uniform distribution between and 1, and the parameter a 
controls the strength of the disorder. We suggest, in analogy with the average length of the optimal 
path, that the distribution of optimal path lengths has a universal form which is controlled by the 
expression — — ^ , where ^oo 

is the optimal path length in strong disorder (a — > oo) and p c is the 
percolation threshold. This relation is supported by numerical simulations for Erdos-Renyi and 
scale-free graphs. We explain this phenomenon by showing explicitly the transition between strong 
disorder and weak disorder at different length scales in a single network. 
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I. INTRODUCTION: 



Many real world systems exhibit a web-like structure and may be treated as "networks" . 
Examples may be found in physics, sociology, biology, and engineering DQD . The function 
of most real world networks is to connect distant nodes, either by transfer of information 
(e.g. the Internet), or through transportation of people and goods (such as networks of 
roads and airlines). In many cases there is a "cost" or a "weight" associated with each link, 
and the larger the weight on a link, the harder it is to traverse this link. In this case, the 
network is called "disordered" or "weighted" 0,0. For example, in the Internet each link 
between two routers has a bandwidth or delay time, in a transportation network some roads 
may have only one lane while others may be highways allowing for large volumes of traffic. 
The average length of the opthna, path ^shortest path", in weighted lattices and 



networks has been extensively studied [4j, 
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In weighted networks it is commonly 



assumed that each link is associated with a weight T{ = exp(arj), where is a random 
number taken from a uniform distribution between and 1, and the parameter a controls 
the strength of the disorder. It has been shown [9] that the length of the optimal path in 
such weighted networks scales as 1(a) ~ N Uopt (where v opt is universal exponent) for small 
system size N, and 1(a) ~ logiV for large systems j^j. More precisely: 



la) ~ 



ap c 



(1) 



where p c is the percolation threshold and ~ N Uopt is the optimal path length for strong 
disorder (a — ► oo). For Erdos-Renyi (ER) graphs v opt = 1/3. For scale-free (SF) networks, 
with a power law degree distribution P(k) ~ k~ x , v opt = (A — 3)/(A — 1) for 3 < A < 4 and 
v opt = 1/3 for A > 4 |4J. The function F(u) is of the form: 

. const if u<l 

F(u) = { . (2) 

I \og(u)/u if u^> I 

In this paper we study the following question: how are the different optimal paths in a 
network distributed? The distribution of the optimal path lengths is especially important in 
communication networks, in which the overall network performance depends on the different 
path lengths between all nodes of the network, and not only the average. A recent work has 



studied the distribution form of shortest path lengths on minimum spanning trees 



10]. 



which 



correspond to optimal paths on networks with large variation in link weights (a — > oo). 



We generalize these results and suggest that the distribution of the optimal path lengths 
has the following scaling form: 

The parameter Z = determines the functional form of the distribution. Relation Q is 
supported by simulations for both ER and SF graphs, including SF graphs with 2 < A < 3, 
for which p c — > with system size N ^j] (Section ITT|) . 

The paper is organized as follows: in Section |H] we show results from simulations for 
various ER and SF graphs. In Section ITTT1 we explain these results and also show that the 
optimal path l op t{a) inside a single network scales differently below and above a characteristic 
length £ = ap c . For I < £ it is like strong disorder, while for I > £ the behavior is like weak 
disorder. 

II. ERDOS-RENYI AND SCALE-FREE GRAPHS: 

We simulate ER graphs with weights on the links for different values of graph size N, 
control parameter a, and average degree (k) (which determines p c = l/(k)) - see Table HI 



We then generate the shortest path tree (SPT) using Dijkstra's algorithm |l2j from some 
randomly chosen root node. Next, we calculate the probability distribution function of the 
shortest (i.e. optimal) path lengths for all nodes in the graph. 

In Fig. [I] we plot £ooP(£, N, a) vs. £/£oo for different values of N, a, and (k). A collapse 
of the curves is seen for all graphs with the same value of Z = — — . 

Figure El shows similar plots for SF graphs - with a degree distribution of the form 

n n 

P(k) ~ k A and with a minimal degree m |2y] [27]. A collapse is obtained for different 
values of N, a, A and m, with A > 3 (see Table ITT]) . 

Next, we study SF networks with 2 < A < 3. In this regime the second moment of the 
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14). For 



degree distribution (k 2 ) diverges, leading to several anomalous properties 

3-A 

example: the percolation threshold approaches zero with system size: p c ~ iV x - 1 — > 0, 
and the optimal path length was found numerically to scale logarithmically (rather than 
polynomially) with iV ^j]. Nevertheless, as can be seen from Fig. Eland Table ITTT1 the optimal 
paths probability distribution for SF networks with 2 < A < 3 exhibits the same collapse 
for different values of iV and a (although its functional form is different than for A > 3). 
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III. DISCUSSION: 

We present evidence that the optimal path is related to percolation 0]. Our present 
numerical results suggest that for a finite disorder parameter a, the optimal path (on average) 
follows the percolation cluster in the network (i.e., links with weight below p c ) up to a typical 
"characteristic length" £ = ap c , before deviating and making a "shortcut" (i.e. crossing a 
link with weight above p c ). For length scales below £ the optimal path behaves as in 
strong disorder and its length is relatively long. The shortcuts have an effect of shortening 
the optimal path length from a polynomial to logarithmic form according to the universal 
function F(u) (Eq. |5J). Thus, the optimal path for finite a can be viewed as consisting 
of "blobs" of size £ in which strong disorder persists. These blobs are interconnected by 
shortcuts, which result in the total path being in weak disorder. 

We next present direct simulations supporting this argument. We calculate the optimal 
path length 1(a) inside a single network, for a given a, and find (Fig. HJ) that it scales 
differently below and above the characteristic length £ = ap c . For each node in the graph 
we find Imini which is the number of links ( "hopcounts" ) alo ng t he shortest path from the 
root to this node without regarding the weight of the link 28[ . In Fig. 0] we plot the 
length of the optimal path 1(a), averaged over all nodes with the same value of l m in for 
different values of a. The figure strongly suggests that 1(a) ~ exp(l min ) for length scales 
below the characteristic length £ = ap c , while for large length scales 1(a) ~ l m in l=M ■ This 
is consistent with our hypothesis that below the characteristic length (£ = ap c ) lmin ~ 

logiV 

and 1(a) ~ iV 1 / 3 , while l m in ~ log N and 1(a) ~ log N above. 

In order to better understand why the distributions of l opt depend on Z according to 
Eq. Q, we suggest the following argument. The optimal path for a —>■ oo, was shown to 
be proportional to iV 1 / 3 for ER graphs and jV^ -3 )/^ -1 ) for SF graphs with 3 < A < 4 Q. 

For finite a the number of shortcuts, or number of blobs, is Z = %^ = — . The deviation of 

£ a Pc 

the optimal path length for finite a from the case of a — >• oo is a function of the number of 
shortcuts. These results explain why the parameter Z = determines the functional form 
of the distribution function of the optimal paths. 
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IV. SUMMARY AND CONCLUSIONS: 



To summarize, we have shown that the optimal path length distribution in weighted 
random graphs has a universal scaling form according to Eq. (j3J). We explain this behavior 
and demonstrate the transition between polynomial to logarithmic behavior of the average 
optimal path in a single graph. Our results are consistent with results found for finite 
dimensional systems (isl . tla . IitI l"!^ : In finite dimension the parameter controlling the 
transition is , where L is the system length and v is the correlation length critical 
exponent (for random graphs v = 1 when calculated in the shortest path metric). This is 
because only the "red bonds" - bonds that if cut would disconnect the percolation cluster [l^ 
- control the transition. 
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TABLE I: Different disordered ER graphs with same value of Z 

Fig.m 



Pc a 



The symbols refer to 



In this method, each node is assigned a number of open "stubs" according to the scale-free 
degree distribution P(k). Then, these stubs are interconnected randomly, thus creating a 
network having the required degree distribution P(k). 
[27] Note that the minimal degree is m = 2 thus ensuring that there exists an infinite cluster for 
any A, and thus < p c < 1. For the case of m = 1 there is almost surely no infinite cluster 



for A > A c ~ 4 (or for a slightly different model, X r = 3.47875 



percolation threshold p c 



(k) 



(fc(fc-i)) 



> 1. See 



rr 



24]! 



resulting in an effective 



23l . |2J] for details. 

28] This is done by using the Breadt h- First-Search (BFS) algorithm jl^ . 
[29] For length scales smaller than £ we have l op t = AN 1 / 3 and l m in = -Bin A, where A and B 
are constants. Thus A = exp (l m in/B) and l op t = Aexp (l m i n /3B). Consequently, we expect 
that: £jp = ^£fe»/3g) = ^ exp ^ lm . n _ sBin^/sB]. We find the best scaling in Fig. H for 
B 
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FIG. 1: (Color online) Optimal path lengths distribution, P(l), for ER networks with (a,b) 
Z = — ^ = 10 and (c,d) Z = 3. (a) and (c) represent the unsealed distributions for Z = 10 
and Z = 3 respectively, while (b) and (d) are the scaled distribution. Different symbols represent 
networks with different characteristics such as size N (which determines ~ average 
degree (k) (which determines p c = l/(k)), and disorder strength a - see Table[l]for details. Results 
were averaged over 1500 realizations. 
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TABLE II: Different disordered SF graphs with same value of Z = ^ „ • The percolation threshold 
was calculated according to: p c = ..fL, . The symbols refer to Fig. [21 
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FIG. 2: (Color online) Optimal path lengths distribution, P(l), for SF networks with (a,b) Z = 
^ „ = 10 and (c,d) Z = 2. (a) and (c) represent the unsealed distributions for Z = 10 and Z = 2 
respectively, while (b) and (d) are the scaled distribution. Different symbols represent networks 
with different characteristics such as size N (which determines ~ N Uopt ), A and m (which 
determine p c ), and disorder strength a - see Table UTI Results were averaged over 250 realizations. 




FIG. 3: (Color online) Optimal path lenghts distribution function for SF graphs with A = 2.5, 
and with Z = — ^ = 10. (a) represents the unsealed distribution for Z = 10 while (b) shows the 
scaled distribution. Different symbols represent graphs with different characteristics such as size iV 
(which determines ~ log(iV) and p c ~ A r ~ 1 / 3 ), and disorder strength a - see Table UTI Results 
were averaged over 1500 realizations. 
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TABLE III: Different disordered SF graphs with A = 2.5 and with same value of Z = — — . Notice 

of p c a 

that p c ~ A -1 / 3 — > for A — > oo. The symbols refer to Fig. |21 
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FIG. 4: (Color online) Transition between different scaling regimes for the optimal path length 
1(a) inside an ER graph with N = 128,000 nodes and (k) = 10. (a) shows the unsealed and 
(b) shows the scaled length of the optimal path 1(a) averaged over all nodes with same value of 
l m in- Different symbols represent different values of the disorder strength a. Fig. (b) shows that 
for length scales £(a) smaller than the "characteristic length", £ = ap c , 1(a) grows exponentially 
relative to the shortest hopcount path l m in (see solid line). This is consistent with 1(a) ~ A 1 / 3 
and l m in ~ log A inside the range of size £ = ap c . For length scales above £ both quantities scale 
as log A. 
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